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ABSTRACT 

Reducing the power on small scales relative to the 'standard' ACDM model alleviates a number of 
possible discrepancies with observations, and is favored by the recent analysis of WMAP plus galaxy 
and Lyman-a forest data. Here, we investigate the epoch of reionization in several models normalized to 
WMAP on large scales, and with sufficiently reduced power on small scales to solve the halo concentration 
and substructure problems. These include a tilted model, the WMAP running-index model, and a 
Warm Dark Matter model. We assume that the Universe was reionized by stellar sources comprised 
of a combination of supermassive (~ 200M Q ) Pop III stars and Pop II stars with a 'normal' IMF. We 
find that in all of these models, structure formation and hence reionization occurs late, certainly at 
redshifts below ten, and more probably at z < 6. This is inconsistent (at 2a) with the determination 
of Zreion — 17 from the WMAP temperature-polarization data and is only marginally consistent with 
SDSS quasar observations. The tension between the galactic-scale observations, which favor low-power 
models, and the early reionization favored by WMAP can only be resolved if the efficiency of Pop III 
star formation is dramatically higher than any current estimate, or if there is an exotic population of 
ionizing sources such as mini-quasars. Otherwise, we may have to live with the standard ACDM power 
spectrum, and solve the small-scale problems in some other way. 

Subject headings: cosmology: theory — galaxies: evolution — intergalactic medium 



1. INTRODUCTION 

The first-year WMAP satellite data have yielded a re- 
markable confirmation of what is coming to be regarded 
as the 'Standard Model' of cosmology: a Hubble parame- 
ter Hq ~ 70 km s _1 Mpc , matter and vacuum densities 
consistent with a flat geometry f2 m ~ 1 — Q\ ~ 0.3, and a 
nearly scale- invariant primordial power spectrum, n s ~ 1, 
at least on large scales. There were some surprises, as 
we will discuss below, but in many ways the WMAP re- 
sults highlight the marked success of standard ACDM in 
reconciling diverse observations from the CMB, to galaxy 
clusters, supernovae, and gravitational lensing. This re- 
markable agreement with data on large scales, however, 
brings into sharper focus several nagging discrepancies on 
small-scales. Specifically, observed galaxy central densi- 
ties appear to be significantly lower than the standard 
ACDM model predicts (Flores & Primack 1994; Moore 
1994; Alam et al. 2002; Swaters et al. 2003; McGaugh et al. 
2003; van den Bosch et al. 2003) and ACDM may also 
produce too much substructure compared with observed 
numbers of satellite galaxies in the Local Group (Klypin 
et al. 1999; Moore et al. 1999). A natural way to relieve 
these problems, while maintaining the large-scale success 
of the model, is to reduce the power on small scales rela- 
tive to large, either by introducing a 'tilt' in the primordial 
power spectrum (n s < 1), as expected in some variants of 
inflation (Alam et al. 2002; Zentner & Bullock 2002), or 
by resorting to Warm Dark Matter (WDM) (Avila- Reese 
et al. 2001; Bode et al. 2001; Alam et al. 2002). In these 
models, collapse happens later, when the Universe was less 
dense, and galaxy halos are naturally less centrally concen- 
trated. Later collapse also alleviates the dwarf problem 

1 Space Telescope Science Institute, Baltimore MD 21218; somerville@stsci.edu; mlivio@stsci.edu 

2 Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138; jbullock@cfa.harvard.edu 

3 Hubble Fellow 



(Colin et al. 2000; Knebe et al. 2002; Bullock & Zentner 
2002; Zentner & Bullock 2003), and the 'angular momen- 
tum catastrophe' — the ongoing problem with producing 
disk galaxies with sufficient specific angular momentum in 
hydrodynamic simulations within the ACDM framework 
(Sommer-Larsen & Dolgov 2001). Interestingly, one of the 
surprises of the WMAP analysis is that it favors (at ~ 2a) 
a model in which the spectral index varies strongly as a 
function of wavenumber, dn/dln/c ~ —0.03 (Spergel et al. 
2003), in precisely the manner needed to cure many of the 
problems on small-scales (Zentner & Bullock 2003). 

The second (and perhaps major) surprise from the 
WMAP report was the detection of a large amplitude sig- 
nal in the TE maps, indicating a large optical depth to 
Thomson scattering (r = 0.17 ± 0.04). The straightfor- 
ward interpretation of this result is that the Universe be- 
came reionized at z rc ion = 17 ± 5 (Kogut et al. 2003), 
rather earlier than many had expected. Several workers 
have recently demonstrated that even within the stan- 
dard ACDM framework, models of reionization require 
rather extreme assumptions in order to produce enough 
early star or quasar formation to reionize the Universe by 
z ~ 17 (Wyithe & Loeb 2003; Haiman & Holder 2003; 
Ciardi et al. 2003a; Sokasian et al. 2003; Cen 2003). Put 
in the context of the small-scale crises facing the standard 
model, and indeed, in view of the running-index model 
favored by the WMAP analysis, we arc left with a para- 
dox. The galactic-scale data favor low-power models pre- 
cisely because they produce late structure formation, but 
the high optical-depth measurement seems instead to fa- 
vor early collapse. In fact, the late reionization implied 
by low-power models would more easily accommodate the 
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simplest interpretation of the Gunn-Peterson troughs ob- 
served in several SDSS quasars at z > 6 (Fan et al. 2001; 
Becker et al. 2001), and the high temperature of the IGM 
at z ~ 4 (Hui & Haiman 2003). 

In this Letter we set out to quantify these concerns us- 
ing a simple yet conservative model that is informed and 
motivated by numerical simulations. Our aim is to com- 
pute conservative upper limits on z ro j on for several mod- 
els that have small-scale power reduced to plausibly solve 
the galactic-scale difficulties, including the running spec- 
tral index model favored by WMAP. In what follows, we 
assume that reionization is due to stellar sources alone, 
as AGN are likely to be unimportant to reionization at 
high rcdshift unless there are exotic objects (such as 'mini- 
quasars') that is completely disjoint from the known pop- 
ulation (Schirber & Bullock 2003). 

2. MODELS AND NORMALIZATION 

In order to focus on the ramifications of reducing the 
small-scale power, we will fix the cosmological parameters 
in all of our models to canonical values: Q m = 1.0 — Oa = 
0.3, fl b h 2 = 0.02, and H = 70 km s" 1 Mpc" 1 . Ad- 
ditional parameters of each of the models are described 
in Table 1. In all cases, we assume that the primor- 
dial power spectrum takes the form P(k) cx k Us , with 
a spectral index that can vary as function of scale as: 
n s = n s (ko) + 0.5dn/dlnfc \n(k/ko) (e.g. Kosowsky & 
Turner 1995; Hannestad et al. 2002). In the second column 
of Table 1, we list n s , evaluated at kg = kwMAP = 0.05 
Mpc -1 and in column 3, we give dn/dln/c. In order to 
make contact with previous work, our specific choice for 
normalizing standard ACDM (with n s — 1.0) is a% = 0.9, 
where as (column 4) is the linear, rms fluctuation am- 
plitude of the power spectrum within spheres of radius 
8h _1 Mpc, evaluated at z = 0. This spectral index and 
normalization are very close to the best-fit parameters de- 
rived from the WMAP data alone, without combination 
with other data sets. In order to ensure consistency on 
large scales, we have chosen to fix the normalization of the 
rest of our models to match that of standard ACDM at 

k = kwMAP- 

The first of our models with reduced small-scale fluctu- 
ations, "TILT", has a s = 0.75 and modest tilt, n s = 0.95, 
as expected in many simple models of inflation (see, e.g. 
Kinney 2003). The specific parameter choice is motivated 
by the work of Zcntner & Bullock (2002; 2003) who con- 
cluded that this normalization/tilt is favored by galaxy ro- 
tation curve data, and may also alleviate the dwarf satellite 
problem without the need for differential feedback. Inter- 
estingly, this parameter choice is also nearly identical to 
the best-fit power-law model in the WMAP joint analysis 
with other CMB data, the 2dF Galaxy Redshift Survey, 
and Lyman a forest data (Spergel et al. 2003). The next 
model, "RSI" , is motivated by the same joint WMAP anal- 
ysis, but represents the case in which they have allowed a 
Running Spectral Index. Their best fit has n s — 0.93 and 
dn/dln/c = —0.03. Our normalization for RSI, as = 0.81, 
is slightly below the quoted WMAP value (by ~ 3.5%) be- 
cause of the simple way we have elected to normalize on 
large scales, but this difference is small compared to the 
quoted uncertainty. This RSI case also does quite well in 
explaining the small-scale observations (Zentner & Bullock 



2003). 

Our final model is computed in the context of a WDM 
scenario, in which the primordial power spectrum is scale- 
invariant, but free-streaming damps fluctuations below a 
characteristic scale set by the WDM particle mass m w : 

R f = 0.2 (0 W dm1i 2 ) 1/3 (^) ~ 4/3 Mpc. (1) 

We calculate the WDM spectra assuming the same flat 
cosmology with fi w = Cl m = 0.3, and use the approximate 
WDM transfer function given by Bardeen et al. (1986): 

P(k) = exp{-kR f - (kR f ) 2 ]P CDM (k). (2) 

We choose m w = 1.5 keV, or Rf ~ 0.027 Mpc, in order 
to suppress power below a mass scale of ~ 10 10 M Q . This 
choice of WDM mass (with n s — 1.0) significantly allevi- 
ates the usual small scale problems (Eke et al. 2001; Alam 
et al. 2002; Zentner & Bullock 2003), yet is not ruled out 
by Ly-a forest data (Narayanan et al. 2000) or by pre- 
WMAP constraints from reionization and early structure 
formation (Barkana et al. 2001). We neglect the finite par- 
ticle velocity dispersion, as well as associated phase-space 
restrictions, which are expected to have negligible conse- 
quence for a 1.5 keV model (Alam et al. 2002; Zentner 
2003). 

3. HALO COLLAPSE RATES 

We show in Fig. 1 the fraction of the total mass in 
the Universe that is in collapsed halos in two different 
mass ranges, computed using the Press-Schechter formal- 
ism (Press & Schechter 1974). In the top panel, we 
show the fraction of mass in halos with mass greater 
than 1.0 x lO 6 h _1 M , but with temperature less than 
10 4 K (conversions between mass and temperature are as 
in Somervillc & Primack 1999). As discussed below, we 
expect gas in these halos (sometimes called "mini-halos" ) 
to be able to cool only via H 2 , and based on the results 
of hydrodynamic simulations by Yoshida et al. (2003), 
they may potentially harbor massive Pop III stars. In 
the bottom panel, we show the mass fraction in halos with 
T v ; r > 10 4 K, which should be able to cool via Hi and form 
Pop II stars with a normal Salpeter-like IMF. We see that 
all models with reduced small scale power have dramati- 
cally decreased mass fractions relative to Standard ACDM 
in halos of the scale expected to be eligible for star forma- 
tion at high redshift. The TILT model shows the smallest 
decrease on all scales, while the WDM model shows the 
most dramatic decrease. 

4. STAR FORMATION AND IONIZING PHOTONS 

A basic requirement for star formation is the presence of 
cold, dense gas. The two primary coolants in the early Uni- 
verse, before the production of heavy elements, are molec- 
ular and atomic hydrogen. Atomic hydrogen is inefficient 
at temperatures below ~ 10 4 K, implying that molecular 
hydrogen was probably the main coolant in the first ha- 
los to collapse at z ~ 20-30, with temperatures of a few 
hundred Kelvin. Recent theoretical work suggests that the 
first stars to form in these primordial halos may have been 
extremely massive, ~ 100-600 M©, but that the efficiency 
was rather low, with less than ~ 1 percent of the available 
gas converted to stars or star clusters (Abel et al. 2000, 
2002; Bromm et al. 2002). 
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The overall efficiency of early star formation is regu- 
lated in mini-halos by a complex interplay between the 
destruction of H 2 by UV photons and its catalysis by X- 
rays (Machacek et al. 2001, 2003; Ricotti et al. 2001, 2002; 
Cen 2002), and in larger halos by photo-evaporation and 
supernova feedback (Ciardi et al. 2000; Somerville 2002; 
Benson et al. 2002). The efficiency of production and mix- 
ing of heavy elements also determines the epoch at which 
the IGM becomes sufficiently polluted with metals to allow 
cooling to lower temperatures and fragmentation, leading 
to a shift from the formation of solely supermassive stars 
to a more normal Salpeter-like IMF (Bromm et al. 2001a). 
For a more detailed discussion of this scenario, see the 
excellent review by Loeb & Barkana (2001), as well as re- 
cent papers by Cen (2002), Haiman & Holder (2003), and 
Wyithc & Locb (2002, 2003). 

We model the global star formation rate density (SFRD) 
by assuming that it is proportional to the rate at which 
gas collapses into halos in a given mass range: 

dF h 

p\ =e*p 6 — (M>M crit ), (3) 

where ^-(M > M cr ;t) is the time derivative of the 
fraction of the total mass in collapsed halos with masses 
greater than M cr i t , obtained from the halo mass function 
-^jj(M,z) given by the Press-Schechter model, and pb is 
the mean density of baryons. We adopt M a - lt = M(T v i r = 
10 4 K) and ej, 1 = 0.1 for 'Pop IF halos, and e™ = 0.002 and 
MlH = 1.0 x K^h^M© for 'Pop IIP halos. In Somerville 
& Livio (2003, SL03), we found that this simple prescrip- 
tion, with these parameter values, agrees well in the red- 
shift range 3 < z < 30 with more detailed semi-analytic 
models and hydrodynamic simulations of Pop II star for- 
mation including photoionization and supernova feedback 
(e.g. Somerville et al. 2001; Springel & Hernquist 2003), 
and also with the general analytic arguments outlined by 
Hernquist & Springel (2002). The global star formation 
rate predicted by this approach is also consistent with ob- 
servational constraints at 3 < z < 6 (SL03). Similarly, 
our Pop III SFRD is in good agreement with the detailed 
numerical hydrodynamic simulations of Pop III formation 
by Yoshida et al. (2003), as also shown explicitly in SL03. 
In both cases, the analytic recipe used here slightly over- 
estimates the star formation rate relative to the more re- 
alistic simulations, thus leading to more optimistic results 
for early reionization. As we do not know when the tran- 
sition from solely supermassive Pop III stars to Pop II 
star formation will occur, we conservatively shut off the 
Pop III mode at z < 6 (when the IGM is known to be 
significantly polluted with metals) . Earlier shut-off times 
will only lead to later reionization, making our conclusions 
even stronger. 

For Pop II stars, we use the results of Leitherer et al. 
(1999) for the number of A < 912 A photons produced by 
low metallicity stars with a "bottom-light" Salpeter IMF. 
At ages less than about 3 million years, these stars pro- 
duce about dn^/dt = 8.9 x 10 46 photons s" 1 Mq 1 . For Pop 
III, we assume that each star produces 1.6 x 10 48 photons 
s^ikfg 1 for a lifetime of 3 million years (Bromm et al. 
2001b). We note that, in the spirit of optimism that we 
have adopted throughout, these values are at the high end 
of the estimates of ionizing photon production and lifetime 



for both populations. 

The cumulative number of ionizing photons per hydro- 
gen atom in the Universe is shown in Fig. 2, for each 
of the four models summarized in Table 1. The ionized 
fraction x e is expected to scale as this quantity times 
/csc /ion/C c i um p (e.g., Spergel et al. 2003), where / ion is the 
number of ionizations per UV photon, / csc is the fraction of 
ionizing photons that escape from the galaxy, and C c i ump is 
the clumping factor, reflecting the dumpiness of the IGM. 
Many analytic and numerical studies have addressed the 
issue of the precise number of ionizing photons per atom 
needed to attain 'overlap', or reionization that is ~ 99 per- 
cent complete (Sokasian et al. 2003; Ciardi et al. 2003b; 
Razoumov et al. 2002; Haiman et al. 2001; Miralda-Escude 
et al. 2000; Gncdin 2000; Madau et al. 1999). This num- 
ber is a function of the local IGM density (thus, indirectly, 
of the redshift of reionization), and of the masses of the 
halos harboring the sources, as more massive sources are 
expected to be in denser environments with higher clump- 
ing factors. For reference, Sokasian et al. (2003) find (in 
the highest resolution numerical simulation of reionization 
by stellar sources to date) that a gross budget (i.e., before 
accounting for the escape fraction) of about 5-20 photons 
per H atom was required to achieve a volume- weighted ion- 
ization fraction of 99%. The clumping factors in numerical 
simulations are known to be potentially overestimated due 
to the limited numerical resolution (Haiman et al. 2001; 
Sokasian et al. 2003), and no numerical simulation to date 
has included the contribution from Pop III stars, which 
are likely to be less clustered. We therefore entertain the 
possibility that as few as ~ 2 ionizing photons per H atom 
may be able to do the job. The redshift at which a total of 
2 or 10 ionizing photons per H atom have been produced 
in each of the models is shown in Fig. 2, and recorded in 
Table 1. We regard these values as bracketing a generous 
but plausible range for the expected redshift of overlap in 
these models. 

One can quickly see that it is difficult to achieve reion- 
ization earlier than z ~ 10 in any of the models with re- 
duced small scale power. In the RSI model favored by 
Spergel et al. (2003), the number of ionizing photons per 
hydrogen atom n 7 reaches the most optimistic range for 
expected overlap (n 7 > 2) at z ~ 8 — only slightly earlier 
than the redshift at which the WDM model reaches this 
range (z ~ 7). Values of n 1 > 10-20, corresponding to 
more typically favored values of the photon escape fraction 
(/esc ~ 0.2), recombination rate, and clumping factor, are 
attained only at z ~ 5.8-4.8 in the models with reduced 
small scale power — thus these models (particularly WDM 
and RSI) may not even be able to reionize the Universe 
early enough for consistency with the SDSS quasar obser- 
vations. As noted above, we have allowed massive Pop III 
star formation to continue until z — 6, although we actu- 
ally expect it to shut off much earlier. This would lead to 
even later reionization, as shown by the short-dashed lines 
in Fig. 2, which show the contribution from Pop II stars 
only. 

4.1. Extreme Parameter Values 

Certainly it may be argued that the IMF and formation 
efficiency of early generations of stars are highly uncertain. 
We now consider how much we would have to vary the free 
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parameters in our simple model to obtain reionization by 
z ~ 17 in the RSI model favored by the analysis of WMAP 
combined with other data (Spergel et al. 2003). These re- 
sults are summarized in Table 2. In model RSI.xl, we 
increase the efficiency of star formation in 'Pop II' halos 
to 100% (corresponding to all baryons forming stars the 
moment they are within a halo) and the assumed mass 
of Pop III stars to 600 M©, the maximum expected value 
suggested by Omukai & Palla (2003), though this is larger 
than the maximum value advocated by Abel et al. (2002). 
Reionization is then expected between z ~ 10.5-8. Only if 
we assume 100% star formation efficiency in Pop III halos 
as well (model RSI.x2) do we find that reionization could 
plausibly occur by z ~ 16-14, in reasonable agreement 
with the WMAP TE results. 

Alternatively, we can leave the star formation efficiency 
parameters at their fiducial values and vary the number 
of ionizing photons produced by each population. For 
example, the IMF of early Pop II stars might be top- 
heavy, because of the higher temperatures and pressures 
and lower metallicities at early times (Larson 1998). More- 
over, Tumlinson et al. (2003) find that zero-metallicity Pop 
II stars with a Salpeter IMF produce only about 50 % 
more hydrogen-ionizing photons than the low-metallicity 
Leitherer et al. (1999) models used here (see also Schaerer 
2003). We find that increasing the number of ionizing 
photons per solar mass of Pop II stars by a factor of 2-10 
can only shift reionization to at best z ro j on ~ 10 (models 
RSI.x3 and RSI.x4). Even with the extreme assumption 
that Pop II stars produce as many ionizing photons per 
unit mass as supermassive Pop III stars (about 20 times 
the fiducial value), we find it is unlikely that the Universe 
could be reionized earlier than z ~ 11 in the context of 
the WMAP RSI model. 

5. DISCUSSION AND CONCLUSIONS 

We have investigated early structure formation in the 
Standard ACDM model and in three models with reduced 
small scale power, each of which alleviates conflicts with 
observations on sub-galactic scales at low redshift. All the 
models we considered are consistent with the WMAP data 
on large scales, and two of them (TILT and RSI) are fa- 
vored by the combined analysis of WMAP and other CMB 
data with Lyman-cv forest and 2dFGRS data (Spergel et al. 
2003). We model star formation using a simple recipe that 
has been calibrated against more detailed semi-analytic 
and numerical hydrodynamic simulations of Pop II and 
massive Pop III star formation (Somerville & Livio 2003; 
Springel & Hernquist 2003; Yoshida et al. 2003). Using 
these fiducial parameter values, and assuming that a gross 
production of between 2-10 ionizing photons per H atom 
is needed to attain overlap, we estimate that the Universe 
would become reionized between z ~ 13 and z ~ 9 with 
the Standard ACDM power spectrum, in agreement with 
several other semi-analytic and numerical studies of reion- 
ization (Gnedin 2000; Razoumov et al. 2002; Haiman & 
Holder 2003; Ciardi et al. 2003a; Sokasian et al. 2003). 
Using the same approach, we find that none of the models 
with reduced small scale power can produce at least two 
ionizing photons per atom before z ~ 9.5, which is about 
1.5 a lower than the epoch of reionization favored by Kogut 
et al. (2003) based on the WMAP TE measurement. More 



plausible values of n 7 ~ 10 are not attained until z < 6 
in the models with reduced power on small scales. The 
WDM and RSI models may have difficulty reionizing the 
Universe even by z > 6, as required by observations of 
high redshift quasars (Fan et al. 2001). Of the three mod- 
els we considered, the model with a fixed tilt n s = 0.95 
(TILT) is the best candidate for obtaining a compromise 
between the requirements of early structure formation and 
observed galaxy central densities and sub-structure at low 
redshift. The RSI model favored by Spergel et al. (2003) is 
only marginally better than WDM in terms of the reioniza- 
tion constraints. It is also worth noting that any attempts 
to push reionization to higher redshift by increasing the 
small scale power (e.g., by adopting a 'red tilt' n > 1 
as suggested by Cen 2003) would further exacerbate the 
problems on sub-galactic scales. 

We find that varying the efficiency of star formation or 
the stellar IMF within a reasonable range of values cannot 
significantly change these conclusions. Only if the effi- 
ciency of star formation is pushed to an extreme upper 
limit (100% of baryons in halos with M > 1 x lO 6 h~ 1 M 
turn instantly into stars) can the RSI model plausibly 
rcionize the Universe by z ~ 16-14. This result is in 
agreement with a similar analysis performed by Haiman 
& Holder (2003). If such extreme star formation were al- 
lowed to continue to lower redshift (even to z ~ 4), such 
a model would be in clear conflict with observations of 
galaxies and the IGM. Even stronger constraints will soon 
be obtained from direct observations of z > 6 galaxies with 
HST and the planned JWST telescope. 

We conclude that if we require that the Universe was 
reionized at least by z ~ 12, within <~ la of the WMAP 
result, then, in the context of the current framework, this 
paradox can be resolved only by adopting one or more of 
the following: (1) the efficiency of Pop III star formation 
is much higher than current theory suggests, (2) there is 
an additional population of ionizing sources (such as mini- 
quasars) at high redshift, (3) the power spectrum has some 
higher-order feature that produces an upturn in power at 
masses just below the scale of dwarf galaxies, or (4) we 
must retain the scale-invariant, n = 1 standard ACDM 
power spectrum, and somehow solve the small-scale prob- 
lems in another way. While none of these solutions is 
particularly attractive, we tend to favor number (4) from 
the context of reionization, although this has its own dif- 
ficulties. If reionization indeed occurred early, and the 
small-scale problems are as robust as they appear, then 
there are significant gaps in our theoretical understanding 
of first light and the structure of galaxies on kpc scales. 
The connection between the two certainly warrants more 
examination. 
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Table 1 
Summary of models 



Model 


n(k W MAp) 


dn/dlnfc 


^8 


m v (keV) 


z(n 7 = 2) 


z(n 7 = 10) 


ACDM 


1.0 


0.0 


0.90 


0.0 


13.3 


8.7 


TILT 


0.95 


0.0 


0.75 


0.0 


9.4 


5.8 


RSI 


0.93 


-0.03 


0.81 


0.0 


7.9 


5.0 


WDM 


1.0 


0.0 


0.89 


1.5 


6.9 


4.8 



Table 2 
rsi model variants 



model 


el 11 


ei 11 


N 11 

iV 7 


z(n 1 = 1) 


z(n 7 = 2) 


z{n 1 = 10) 


RSI.xl 


1.0 


0.006 


X 1 


11.4 


10.5 


8.6 


RSI.x2 


1.0 


1.0 


X 1 


16.5 


15.8 


13.9 


RSI.x3 


0.1 


0.002 


x 2 


9.6 


8.7 


6.0 


RSI.x4 


0.1 


0.002 


x 10 


11.14 


10.3 


8.4 


RSI.x5 


0.1 


0.002 




11.7 


11.0 


9.1 



The Epoch of Reionization 
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1 — i — i — r 




5 10 15 20 25 30 

redshift 

Fig. 1. — The fraction of mass in collapsed dark matter halos as a function of redshift for the four power spectra summarized in 
Table 1. The top panel shows the mass fraction in halos with mass greater than 1.0 x lO 6 h _1 M0 but virial temperature T vlI < 10 4 
K. In the absence of metals, gas in these halos can cool primarily by molecular hydrogen, and we associate them with possible sites 
of massive (200M Q ) Pop III star formation. The bottom panel shows the mass fraction in halos with T v ; r > 10 4 K. We associate 
these halos with H/ cooling and Pop II star formation with a 'normal' IMF. The impact of reduced small scale power on early 
star formation on both of these mass scales is dramatic. 
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5 10 15 20 25 30 5 10 15 20 25 30 

redshift 

Fig. 2. — The cumulative number of hydrogen-ionizing photons per hydrogen atom produced, as a function of redshift, for the 
four initial power spectra. Long dashed lines show the contribution from Pop III stars, short dashed lines show the contribution 
from Pop II stars, and solid lines show the total. Straight dotted lines indicate the redshift at which 2 and 10 ionizing photons 
per atom have been produced, and roughly bracket the range of redshifts when the overlap of ionized regions is expected to occur. 
In all three models with reduced small scale power, overlap occurs later than z ~ 10. This may be compared with the estimate of 
Zi-cion = 17 ± 5 from the WMAP temperature-polarization results (Kogut et al. 2003). 



